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Abstract 



In the attempt to fully profile a Z' boson accessible at the barge Hadron Collider 
Q\ (LHC), we study the sensitivity of di-lepton (for the electron, muon and tauon cases) 

and di-quark (for the case of the heavy flavours, t and, possibly, b) samples to the 
nature of the new gauge state, for a one-dimensional class of non-exotic Z' bosons. As- 
suming realistic final state reconstruction efficiencies and error estimates, we find that, 
depending on the CERN collider energy and luminosity, the best chances of extracting 
the Z 1 quantum numbers occur when two or more of these channels are simultaneously 
explored, as none of them separately enables one to fully probe the parameter spaces of 
the aforementioned models. Effects of Standard Model (SM) background as well inter- 
^ ferences between this and the various Z' signals have been accounted for. A complete 

study of cross sections and asymmetries (both spatial and spin ones) makes clear the 
need for complementarity, especially for their disentanglement over the full parameter 
space. 



1 Introduction 

So-called Z' bosons, heavier versions of the SM Z gauge boson, appear in many Beyond the 
SM (BSM) scenarios: e.g., in Little Higgs models p], in frameworks with extra spacial dimen- 
sions [2 J as well as in Grand Unified Theories (GUTs) such as 50(10) [3j (e.g., in Left-Right 
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symmetric models) and Eq jl] types. Furthermore, in some realisations of Supersymmetry 
(SUSY) [5] and in Hidden Valley models [6J, they manifest themselves rather subtly, as 
messengers between the SM and the hidden sectors corresponding to such scenarios. 

Such Z' states are best searched for at hadron colliders through the Drell-Yan process, 
i.e., a di-lepton signature proceeding via pp(p) — >■ (7, Z, Z') — >■ £ + £~, where £ — e, fi. In fact, 
the most stringent limits on Z's at both the Tevatron and LHC emerge from these final states. 
Currently, assuming a sequential Z', it should be noted that the Tevatron places limits on 
Mz>, at around 1 TeV [7j while the LHC sets about 2.3 TeV [8]. Recall that a sequential Z' 
is a state with generic mass and same coupling to the SM particles as the Z boson, so that 
limits in the aforementioned models must be obtained by rescaling the results for such a Z' 
bosor^] While the above DY signature is clean experimentally and theoretical uncertainties 
are well under control, it is by no means the only one accessible at the CERN machine. In 
particular, heavy quark final states, i.e., those accessible via pp(p) — >■ (7, Z, Z') — > QQ, where 
Q — b,t, also yield detectable rates. Further, they can both be used to profile the Z' boson 
(that is, to measure its quantum numbers), alongside the DY channel, at the LHC (see, e.g., 
[U]). The case has been made recently for Q = t by ourselves [TU], building on previous results 
[TT] . Here, one exploits the fact that the top (anti)quark pairs decay before hadronising, so 
that their spin properties are effectively transmitted to their decay products. This enables 
one in turn to access the vector and axial couplings of the Z' from which the pair originates. 
Further, recall that the electromagnetic charge of the top can be tagged (via a lepton and/or 
a 6-jet) [12], thanks to which one can define suitable spatial/spin asymmetries which turn 
out to be particularly effective to pin down said couplings of the new gauge boson. On the 
same footing, r leptons are also short lived so that their spin information is imprinted in 
their decay product kinematics [13] . Similarly to top quarks, r polarisation can be measured 
by means of its one-prong decays and analyses exist at the LHC already [H]. 

The case of Q = b brings further complications. The displaced vertex tagging allows 
for the charge measurement of the 6-quark in the jet, which is needed in the definition of 
the spatial asymmetries in 6-quark final states, which therefore seem feasible quantities to 
deal with. Here the main issue comes from the 6-tagging performances at high which we 



discuss in section [2.4.4| Nonetheless, one might hope the tagging efficiency to improve due to 
upgrades in the micro-vertex detectors at both ATLAS and CMS for the \fs = 14 TeV run 
(see, e.g., [L7]). Regarding spin asymmetries, they certainly can be defined also for 6-quarks. 
If the 6-quark hadronises before decaying, it has also been shown that its spin informations 
can be transmitted to 6-hadrons (e.g., A;,), albeit with some dilution factors [IE]. Then, the 
semi-leptonic decays of the latter would enable spin measurements, as pioneered at LEP [19J. 
It is not clear however if the same measurement is possible at the LHC, due to the higher 
boost of the emerging lepton and the narrower 6-jet cone. Despite being nowadays technically 
challenging and still uncertain for the very near future, as better vertex detectors will be 
available, it is still interesting to study spin asymmetries for 6-quarks in view of when such 
measurements will be available. We will therefore present them in the spirit of suggested 
observables. 



however, notice that this implicitly assumes that the Z 1 cannot decay into additional extra matter 
present in the model spectrum. We will dwell on this aspect later on in the paper. 
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It is the purpose of this paper to compare the yield of £ + £~ (£ = e//i, r), bb and tt final 
states produced at the LHC in presence of a Z' boson and to assess the machine ability to 
profile the latter in both standard kinematic variables as well as spatial/spin asymmetries, by 
continuously scanning the parameter space of a one-dimensional class of Z' models, defined 
in terms of the Z' coupling and mixing parameters, over which we will define benchmark 
points and lines amenable to experimental investigation. In doing so we will borrow some 
(but not all) of the standard Z' benchmarks from, e.g., j20] and [21] . as well as define new 
ones. 

The plan of the paper is as follows. In the next section we describe the model, outline 
our calculation and define the observables to be studied. In sect. [3] we report and comment 
on our results. Sect. [4] finally presents our conclusions and outlook. 

2 Framework 

We describe here the model we studied, fixing conventions and relevant features, the details 
of the code used and the variables that have been analysed. 

2.1 The minimal Z' model 

The general class of models we are going to study is the one defined by the so-called non- 
exotic minimal Z' models [211 E21 1231 12H ESI Following the convention of Ref. [2IJ, we 
summarise here just the relevant parts, i.e., the gauge and neutrino sectors, and we refer the 
reader to the latter publication for a complete description of the model. 

When the global U{1) E -l symmetry of the SM is promoted to local, the new Z' gauge 
boson will mix with the SM Z boson. In the field basis in which the gauge kinetic terms are 
diagonal, the covariant derivative ig^J 

D^ = d, + ig s T a G« + igT a W; + i 9l YB^ + i(gY + g'^B-L^ . (1) 

In our bottom-up approach, we will not require gauge unification at some specific, yet 
arbitrary, energy scale, fixing the conditions for the extra gauge couplings at that scale. 
Therefore, in this model, the gauge couplings g and g[ are free parameters. To better 
understand their meaning, let us focus on eq. ([TJ. The last term of the covariant derivative 
can be re- written defining an effective coupling Y E and an effective charge g^. 

g E Y E = gY + g[Y B _ L . (2) 

As any other parameter in the Lagrangian, g and g[ are running parameters [22j [2U 1251 12?]. 
therefore their values ought to be defined at some scale. A discrete set of popular Z' models 
(see, e.g., Refs. [231 EE]) can be recovered by a suitable definition of both g and g[. 

2 In all generality, Abelian field strengths tend to mix. The kinetic term can be diagonalised with a 
GL(2, R) transformation, leading to the form of the covariant derivative in eq. ([!]), where the mixing between 
the two f/(l) factors becomes once again evident [2211121] . 
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Furthermore, three right-handed neutrinos are required by the anomaly-free conditions, 
naturally implementing a type-I seesaw mechanism via Yukawa interactions with the SM 
Higgs field and the new Higgs field needed for the U(1)b-l symmetry breaking to give the 
Z' boson a mass. A general feature of the seesaw mechanism is that the mass eigenstates, 
called "light" (u{) and "heavy" neutrinos, are Majorana particles. Because of this, 
neutrinos only have axial couplings to the neutral gauge bosons. One of the aim of this 
paper is to investigate the impact of such heavy neutrinos on the observables under study. 
We will therefore restrict ourselves to two opposite scenarios: a "decoupled" case, with 
heavy neutrinos much heavier than (half of the) Z' mass, thereby disallowing (even off-shell) 
Z' — > VhVh decays, and a "very light" case, where heavy neutrinos are much lighter than the 
Z' itself (e.g., m Vh = 50 GeV, compatible with LEP limits). All possible intermediate cases 
will therefore lie somewhere in between. 

2.2 Structure of the chiral couplings 

To study the asymmetries it is important to understand the chiral structure of the couplings 
of the Z' gauge boson to fermions, determined by the co variant derivative of eq. (JIJ. 

Due to the mixing between the Z and Z ! gauge bosons, the couplings of the Z' to fermions 
are a function of g[ and g (notice that also 8', the angle parametrising the Z-Z' mass mixing, 
is function of the gauge couplings). In all generality, such couplings can be separated into 
vector and axial components, Cy(g[,'g) and C^(g[,g), respectively, with an interaction term: 



in which we have re-absorbed the gauge couplings in the definition of the vector and axial 
components. 

Because of universal couplings, we can render explicit the chiral structure of the Z' 
couplings to up-type quarks, down-type quarks and to charged leptons, independently of the 
fermion generation, as well as the partial decay widths in the approximation of a massless 
final state. 

Up-type quarks (u): 




(3) 



C C w (Ag', + 5g) S w + e S' (3 - 8 



V) 



(4) 






(5) 
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Down-type quarks (d): 



nd( , C C w j-4g[ + g)S w + eS' (-3 + 4 S 2 W ) 

^vKdi^g) — 77777 q > \<) 

cm,?) = c,5< r;V e5 ' - < 8 > 

4 L<w >JW 

T d (g[,g) = T(Z> -+ dd) ~ ^ (V, - £) 2 + ^) • (9) 

Charged leptons (£): 

^//^ _ - C 'Cw M + m S w + eS' (1 -4 5%,) 

WlS'ufi'J — 777 o j l iU J 

Ca{9i,9) = 777 q > I 11 ) 

r^S) = r(z' £+r) ~ ^ f 1 (4^ + s^) 2 + l? 2 ) . (12) 

In the previous formulas, S'(C) = sin0'(cos#') and Sw{Ow) is the sine(cosine) of the weak 
mixing angle. Due to the smallness of 9', as required from LEP constraints [21], sin#' <C 1 
and cos 9' ~ 1 have been adopted in the simplified evaluation of the partial widths. 

For the evaluation of the total width, we need to include also the light and heavy neu- 
trinos. Their partial widths read as, respectively, 

r^,3f) = r(^^^) ~ ^\( 2 9i + 9)\ (is) 

T Vh {g[)=T{Z' ^v h v h ) ~ ^(g'lf- (14) 

The extra overall factor 1/2 in the previous formulas reflect the Majorana nature of the 
neutrinos, that have only axial couplings. The total width is then given by the sum of the 
partial contributions 

T z , = 3 (r u + V d + T e + T Ul + Y Uh ) , (15) 

where all fermions have been taken massless in comparison to the O(TeV) mass for the Z' 
bosonE 



2.3 Benchmark models 

From the previous equations, benchmark models can be selected for specific features. Com- 
monly employed in the literature are the 'pure' B — L model, that is defined by the condition 
g = (implying no mixing at the tree level between the B — L Z' , sometimes denoted as 
Z'b-l-i an d the SM Z, also denoted as Zsm)', the U(1)r model, for which Left-Handed (LH) 
fermion charges vanish (recovered here by the condition g = — 2g^); the 'S'O(10)-inspired' 

3 In the case of decoupled heavy neutrinos one of course has T Uh = 0. 
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U(l) x model, that in our notation reads g = — \g'\ (the only orthogonal U(l) extension of 
the SM hypercharge) . 

The smallness of 9' allows us to neglect terms proportional to sin#' in eqs. (|4j) — ( 11), 



which elucidates the qualitative features of these common benchmark models. The x model 
will have C v ~ 0; the R model will have equal couplings (in absolute value) to all the 
fermions (\C V \ ~ \C A \ ~ \C V \ ~ \C A \ ~ \C V \ ~ \C A \) while the pure B — L model (for which 
g = S m9' = 0) will have C\ = C d A = C A = 0. 

Spatial asymmetries have the following features: they vanish at the Z' peak if any of 
the chiral couplings are zero and are very small around it mainly due to the interference 
with the SM background. In our framework, the axial couplings are proportional to g, and 
hence vanish identically only in the pure B — L model. Regarding the vector couplings, 
benchmark scenarios can now be identified by requiring some of them to vanish. Beside the 
U(l) x model, already known in the literature, the ft ("B-not") and the p ("L-not") models 
are here defined^ These three models are characterised by the vanishing vectorial coupling 
to the up-type quarks, down-type quarks and charged leptons, respectively, and are defined 
by the following relations: 

U{l) x model: 9 = ~t 9i ~+ C v ~ , (16) 
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c v 


~ 0, 


— > 


Cy 


~ 0, 


— >■ 


Cy 


~ 0. 



(17) 

L model: g = ~g[ -)• C v ~ . (18) 

The last two models are newly proposed here. These models will have trivial asymmetries 
in the corresponding final state (the one with negligible vector coupling) and non- vanishing 
ones in the other final states. Only the pure B — L model will have trivial asymmetries in 
all final states. Finally, the R model will have almost the same values for the asymmetries 
in all final states. 

On the scenario line characterised by a specific relation between the gauge couplings, we 
select specific values for the latter to recover well known benchmark models, some of which 
correspond to those described in Ref. [20]. It is important to note here that only models 
that are anomaly-free with the SM particle content plus Right-Handed (RH) neutrinos are 
included in our framework. Finally, to recover the S'O(10)-inspired U(l) x point, its gauge 
couplings have been rescaled by the usual factor to account for gauge coupling unification 
(a/3/8 in our notation). 



2.4 Calculation and variables 

The code exploited for our study of the asymmetries is based on helicity amplitudes, defined 
through the HELAS subroutines |30j, and built up by means of MadGraph [31]. Initial 
state quarks have been taken as massless while for the (anti)top state we have taken m t = 
172.9 as pole mass. For the 6-(anti)quark we have taken mj, = 4.95 GeV (again, as pole 

4 Not to be confused with the models in which the baryon number and the lepton number are separately 
gauged. 
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mass). The electron and muons were taken as massless while for the tauon we have used 
m T = 1.77 GeV. The Parton Distribution Functions (PDFs) exploited were CTEQ6L1 [32J, 
with factorisation/renormalisation scale set to Q — fi — Mz>- VEGAS [33] was used for the 
mult i- dimensional numerical integrations. A separate program was also used for part of the 
analysis, based on CalcHEP [31], wherein the model has been independently implemented 
via the LanHEP module [35] and the FeynRules [3B] package. The model files can be found 
on the HEPMDB database |37J and on the FeynRules website [38] . 

For all final states, / = £, Q, a cut in the invariant mass window around the Z' peak was 
performed 

\m(fj) -M z ,\< 100 GeV, (19) 

which enhances the Z' peak with respect to the SM background. Regarding the 6-quark final 
state, this is still not sufficient to isolate the resonance. Consequently, a further selection 
has been implemented: 

Prip) > 300 GeV. (20) 

The charge/spin variables that we are going to study have been described in Ref. [10] 
and we summarise here their salient features. The dependence on the chiral couplings of 
the asymmetries can be expressed analytically, using helicity formulae from Ref. [10] (also 
derived independently with the guidance of [H]), for a neutral gauge boson exchanged in 
the s-channel. 



2.4.1 Charge asymmetry 

Charge or spatial asymmetry is a measure of the symmetry of a process under charge conju- 
gation. Due to the CP invariance of the neutral current interactions, this translates into an 
angular asymmetry at the matrix element level. It can only be generated from the qq initial 
state due to the symmetry of the di-gluon system and it occurs dominantly from NLO QCD 
interferences, as described in detail in [39J in the case of final state quarks, while for leptons 
the contribution comes only from the EW sector. 

The symmetric pp initial state at the LHC necessitates a more suitable definition of 
such an observable compared to the well-known top quark forward-backward asymmetry 
employed at the Tevatron. Several possibilities exist, though as investigated in [10J, the 
spatial asymmetry that delivers the higher sensitivity is the rapidity dependent forward- 
backward asymmetry, Arfb- It uses the rapidity difference of the final state fermion pair, 
Ay = \yf \ — \yf\, and enhances the qq initial state parton luminosity via a cut on the rapidity 
of the fermion-anti-fermion system, yjj, 



A 



RFB 



N(Ay > 0) - N(Ay < 0) 



(21) 

\Vff\>Vff 



oc 



N(Ay > 0) + N(Ay < 0) 

C^yCi&y . (22) 



In this paper, we used y c ^j = 0.5. It is clear that this observable can only be generated by 
a Z' boson if its vector (Cy) and axial {Ca) couplings to both the initial (i) and final (/) 
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state fermions are non vanishing. We will study A^fb for light leptons [I = e, /x) and for 
the heavy quark final states. 



2.4.2 Spin asymmetries 

Spin asymmetries, on the other hand, focus on the helicity structure of the final state fermions 
and, when such properties are measurable, display interesting dependences on the chiral 
structure of the Z' boson couplings. The helicity of a final state can only be experimentally 
determined for a decaying final state, where the asymmetries are extracted as coefficients in 
the angular distribution of its decay products. This is described for the case of top quarks 
in Ref. [12], but can be easily extended to other weakly decaying particles. 

We define two such asymmetries. The main observable we consider is the polarisation or 
single spin asymmetry, Al, defined as follows: 

A N(- -) + N(- +) - N(+, +) - iV(+, -) 

a l = ^ , (2d) 

^Total 

cx C{C^ ([C^yf + {C\f) , (24) 

where (5 = ^Jl — 4 m'j/s and iV denotes the number of observed events and its first (second) 

argument corresponds to the helicity of the final state particle(anti-particle). It singles out 
one final state particle, comparing the number of positive and negative helicities, while 
summing over the helicity of the other anti-particle (or vice versa). This observable is 
proportional to the product of the vector and axial couplings of the final state only and is 
therefore additionally sensitive to their relative sign, a unique feature among asymmetries 
and cross section observables. This is equivalent to a sensitivity to the relative right- or 
left-handedness of the Z' couplings. As described in the introduction, it can be measured 
now for t and t-quark final states while for 6-quarks an improved vertex detector would 
probably be needed, so that we only present the latter case here as a suggested analysis. 

In the case of appreciably massive final states, i.e., the top quark only, the spin correlation 
or double spin asymmetry is accessible. This observable relies on the helicity flipping of 
one of the final state particles, whose amplitude is proportional to mf/yf§, where y/S is 
the partonic centre-of-mass (CM) energy, and gives the proportion of like-sign final states 
against the opposite-sign states, 

A iV(+, +) + N(- -) - N(+, -) - N(- +) 

All - 77 , (/oj 

^ Total 

oc (3 (Cl) 2 /3 2 + (C<£) 2 (2 + /3 2 )) ((Q) 2 + (Q) 2 ) . (26) 

This observable depends only on the square of the couplings in a similar way to the total 
cross section. In the massless limit, — > 1 and All becomes maximal, making it a relevant 
quantity to measure only in the tt final state. 
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2.4.3 Asymmetries and Z' couplings 

The ultimate scope of profiling a Z 1 boson is that of measuring its couplings to fermions. 
Although it is beyond the scope of this paper, we still would like to make some comments. 
In the minimal case considered here, 5 independent parameters exist (q^, u R) d R , £ L , e R , 
where here q^ and £l identify the quark and lepton doublets, respectively), so that 5 in- 
dependent measurements are required. In Ref. [15], it is shown that a set of 4 coefficients 
(c q , e q ; q = u, d) that are functions of the Z' couplings can be extracted from 4 observables 
related to charge asymmetry and total cross section in the light lepton channel (electrons 
and muons), by restricting the kinematical domain of the differential cross section. It is 
observed that a degeneracy exists between the dependence on leptonic and quark couplings 
in this analysis, stemming from the minimal assumption of 5 independent parameters and 
only 4 independent observables. Including the total Z 1 width does not seem to fix the prob- 
lem, since it enlarges the set of independent parameters to include all the fermions, such as 
the SM light neutrinos and all the coupled exotic particles, in our case the heavy neutrinos. 
Instead, spin observables, being based on polarised amplitudes, cannot be obtained from the 
fully differential cross section and hence can be considered. Staying on the leptonic final 
state, the tau polarisation is the first case to look at. It has a different dependence on the 



quark and lepton couplings than c q , e q . Recasting eq. (23), it can be written in the spirit 
of [15] as a new coefficient, 

/ ? =„ >d cx (q 2 R + q 2 L )(e 2 R -el), (27) 

making manifest the aforementioned sensitivity to the 'handedness' of the Z' couplings. 
These coefficients can be extracted from tau polarisation measurements in a similar way 
as the c q and e q coefficients are from the differential cross section, i.e., by splitting the 
kinematical domain of A T L in 2 independent regions to then be fitted. It is evident that 
the same quark-lepton degeneracy exists and in fact f q=u ,d turns out not to be linearly 
independent from the c and e coefficients when ul = d^. 

f q =u,d oc c q 6u _ 6d . (28) 

Although the tau polarisation does not allow for the extraction of the Z' couplings within 
the 5 parameter minimal assumption, it still proves to be a useful quantity to measure. 
First of all, it can be used as a test of universality, having better sensitivity than the charge 
asymmetries in a tau final state. Secondly, it provides an extra set of constraints on the 
other coefficients in order to improve the quality of a fit of the couplings. The ratios 

fu _ fd _ e R ~ e L (29) 

should be equal and are independent of the quark couplings and could serve to reduce some 
systematics, such as the one originating from the PDFs. To break the degeneracy allowing for 
a direct fit of the couplings, one must then consider observables in non-leptonic final states, 
e.g., top polarisation or charge asymmetry. For the longer term, 6-quark final states may turn 
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out to be useful too. Moving away from the minimal assumption permitting ul 7^ di (but 
still requiring universality), as could occur with Z's arising from more general gauge group 
extensions, restores the linear independence of one of the tau polarisation coefficients. In this 
case one still requires the use of an alternative final state to complete the set of independent 
observables needed to fit directly to the 6 (or more) couplings, to that the exploitation of 
top (and possibly bottom) final states is then called for again. 

2.4.4 Efficiencies and uncertainties 

To be able to quantitatively address the distinguishability among the various models and 
the SM background, we associate a statistical error to each observable. Although in this 
work we only estimate statistical uncertainties, systematics may also be important J42J 143] 
though the mass window selection is expected to milden their actual contribution. However, 
their inclusion would require detailed detector simulations which are beyond the scope of 
this paper. In this respect, it should further be noted that, by the time the LHC reaches the 
14 TeV stage, where our most interesting results are applicable, systematics will be much 
better understood than at present. 

Given that an asymmetry is defined in terms of the number of events measured in some 
generic 'forward' (Np) and 'backward' (Nb) directions (this is also true for spin asymmetries), 
the statistical error is evaluated by propagating the Poisson error on each measured quantity 
(i.e., 8N F ( B ) = y/Np(B))- Per given integrated luminosity £, the measured number of events 
will be N = eCa, a being the cross section, yielding an uncertainty on the asymmetry A of 



Here, e corresponds to the reconstruction efficiency of the final state system. For the tt 
system, we take e t = 10% as in [TO] , considering all possible decay channels, based on 
efficiencies quoted by recent experimental papers [44J as well as estimates from Monte Carlo 
studies in previous works [45j . For the bb system, we take £& = 5%. This number is based on 
an LHC 6-tagging efficiency of 50%, given by summing 1-tag and 2-tags efficiencies at high 
Pt [16] . to extract the 6-quark decay mode from the di-jet sample. The latter reference shows 
that, while the 2-tags do not yield measurable rates, the 1-tag alone seems insufficient, given 
the almost equal rate of mistags. However, it is shown in [16] that additional requirements 
on a muon from the 6-jet is capable of improving both the tagging and the tagging-over- 
mistag efficiencies for 6-quarks up to intermediate jet pt while no information is available 
for very energetic jets. Regarding the 6-hadron helicity, there is no result currently available 
at the LHC, thus we take the values available from LEP [T51 [T5] . quantifiable in an overall 
10% efficiency, as an estimate to quantitatively discuss spin observables. By considering 
this extra efficiency factor, we are in fact also further selecting 6-jets with prompt hard 
leptons. We base our overall reconstruction efficiency estimate on all of these considerations 
hoping that, for the time scales at which this paper could be relevant, the situation will have 
improved as suggested by the proposed new generation of micro-vertex detector available 
for the y/s = 14 TeV run, see, e.g., Ref [17J. However, we remind the reader that these 




(30) 
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observables are just a suggested measure and might not be available at the LHC, although 
the proposed upgrades for micro-vertex detectors may well help to extract them. For the 
system, we distinguish between light leptons (e, /x), for which e e ^ = 90% [4"7| . and the 
taus, for which e T = 10%. This latter number is based on a recent ATLAS measurement in 
1 -prong decays (25% of tau's Branching Ratios (BRs)), its reconstruction efficiency of 60% 
and a 20% loss of sensitivity due to detector effects [14]. We remark here again that light 
leptons will be employed for charge asymmetries only, whilst we will make use of taus to 
define the spin asymmetry Al in the leptonic final state. 

All figures that will be shown in the following are for the LHC at \fs = 14 TeV assuming 
100 fb _1 of integrated luminosity, for a Z' mass of 2.5 TeV and decoupled heavy neutrinos. 
In some cases, observables are evaluated also for 'light' heavy neutrinos of 50 GeV mass to 
estimate their impact. Values of an observable in the (g[,g) coupling plane can be compared 
to SM predictions corresponding to the point (0, 0). For the signal, the uncertainty for that 
value is considered to confirm visibility over the background in that channel. Also shown are 
surface plots of the significance of the observable with respect to the SM prediction, defined 

s = W)-^)l (31) 
y/6A(l) 2 + 6A(2) 2 ' 

where A(l) and A(2) denote the prediction for an observable by two different hypotheses 
and 5 refers to the statistical uncertainty. In this case, the statistical uncertainty of the 
SM prediction is taken to be zero because its very low cross section at such high invariant 
masses (especially true for the leptonic final state) would artificially dilute the significance of 
asymmetries in that region. Therefore, in the hypothesis of observing a Z' boson, we prefer 
to consider the SM prediction as a reference value only. 



3 Results 

We start by presenting the most up-to-date exclusions from direct searches, from the CMS 
analysis at ~ 5 fb _1 [H] , followed by a description of the features pertaining to each possible 
final state, i.e., BRs, cross sections at the LHC for yfs = 14 TeV (as the scope of 7 and 8 
TeV stages of the CERN machine is rather minimal for our purposes) and total Z' width. 
Notice that ATLAS has also published an analysis for ~ 5 fb _1 [8], but their limits are less 
tight than the CMS ones. Therefore, we will not present them here. We then move on to 
profiling the Z' boson analysing its BRs, cross section and asymmetries in given final states. 



3.1 Exclusion limits 

To start with, we present (see figure [T| the most recent 95% (Confidence Level) CL exclusions 
at the LHC in the {g[,g) plane (first in [IE]), based on the CMS data at y/s = 7 TeV for the 
combination of 4.7(9) fb _1 in the electron (muon) channels [8]. 

Table [I] collects the maximum allowed g[ coupling per given Z' boson mass for the various 
benchmark model of interest (see section 2.3). 
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Figure 1: Z' exclusions from CMS latest data, at a/s = 7 TeV for the combination of 4.7 
fb _1 in the electron channel and of 4.9 fb _1 in the muon channel. The dotted black lines 



refers to the main benchmark models of this analysis (see section 2.3). 
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Table 1: Maximum g[ allowed at 95% CL for the benchmark lines of section 2.3 



3.2 BRs, total widths and production cross sections 

We start by analysing the BRs of the individual final states: up- and down-type quarks, 
charged leptons and light and heavy neutrinos. Minima and maxima are easily found by 
minimising the respective BR. We chose to evaluate them as a function of g: 



d 



(32) 



so to obtain relations of the type g = g(g[). The results with maximum and minimum 
expected BRs are collected in table [2] and are displayed in figure [2] 

We recognise some features: the U(1)r scenario has the lowest BR into light neutrinos, 
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Table 2: g = g(g[) that delivers the maximum and minimum BRs and their extreme values 
(per each generation), for (top ) decoupled and (bottom) very light heavy neutrinos. 



regardless of the heavy neutrino masses. This is simply because light neutrinos are predom- 
inantly LH fermions, hence they do not couple to the Z' boson in this scenario. When the 
heavy neutrinos are light, the maximum of the BR(Z' —y V\Vi) is in the pure B — L model. 
Regarding the heavy neutrinos, they are predominantly RH fields (as already mentioned), 
hence with zero hypercharge. Consequently, their partial width is independent of g, as in 
eq. (14): on the one side, the minimum BR is for g[ —y (when they decouple), on the other 
side the maximum BR is for the minimum total Z' width. By direct calculation, we find that 
the total width of eq. (15) is minimised in the U(l) x scenario (figure [3]), as also indicated in 
table [2] (corresponding to the maximum BR into heavy neutrinos). 

Figure |3] (left) shows the production cross sections for two different values of the Z' mass. 
For large values of the couplings, the production cross section can be up to 0(1 pb), in 
regions of parameter space not excluded experimentally, in accordance with figure [TJ For 
very small albeit non-vanishing couplings, the production cross section for a Z' boson is 
easily bigger than tens of fb. These values are some orders of magnitude above the expected 
SM background when leptons in the final state are considered, as we will see in the following 
section. 

Heavy neutrinos can have a strong impact on the Z' features, altering the lineshape in 
a measurable way. First of all, the total width is augmented by up to a factor 30%, when 
comparing same Z' boson mass and same values for the gauge couplings. A similar factor is 
also the maximum BR of the Z 1 boson into heavy neutrinos (once summed over generations). 
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Figure 2: BRs per each generation (top row) for Z' — > e + e~ and for Z' — > bb, (middle row) 
for Z' — > ti and for Z' — > uiUi, and (bottom row) for Z' — > u^Uh, in the (g[, g) plane for two 
extreme neutrino masses, i.e., for decoupled neutrinos and for very light neutrinos {m Uh = 50 
GeV). The benchmark models are highlighted as (red) dots on the (dotted) scenario lines. 
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Figure 3: Total production cross section (for y/s = 14 TeV) (left) and total width (right) for 
Mz> = 2.5 TeV in the (g[, g) plane for two extreme neutrino masses, i.e., for decoupled neu- 
trinos and for very light neutrinos {m Vh = 50 GeV). The benchmark models are highlighted 
as (red) dots on the (dotted) scenario lines. 

3.3 Total event rates 

It is instructive to present at this point the overall event rates in each channel. When 
including the SM background, the total process reads 

p,p^ 7 , Z, Z'^fJ. (33) 

have been applied 



To evaluate the total cross sections, the cuts described in eqs.(19)- 
to enhance the signal. The numerical results are depicted in figure [4 These are essentially 
convolutions of the BRs and production cross section detailed in figures [2] and [3] and including 
small interference effects with the SM background. These figures highlight the different areas 
of parameter space favoured by each final state and relate directly to the magnitude of the 
statistical uncertainties in its asymmetries (to be studied below). 

As previously observed, the total rate when leptons are considered is above the fb level 
for most of the parameter space. The SM contribution is 0.032 fb, some orders of magnitude 
below the signal also when light heavy neutrinos are considered. For the heavy quarks, 
instead, the SM background is of the same order of the signal, as the former is largely due 
to QCD while the latter is (despite being resonant) an EW process. We further observe 
that the total event rate at the LHC for \fs = 14 TeV can be modified when the decay into 
heavy neutrinos is allowed. In particular, it can diminish by up to 40%, 30%, and 20% when 
considering pp — » , pp — > bb, and pp —> it, respectively. 

If the presence of light heavy neutrinos lowers the rate for a particular final state, directly 
increasing its relative error and therefore the error of the asymmetries, we expect this to 
directly influence the central value of the asymmetries too. Given that the SM background 
is not altered by the presence of heavy neutrinos, the net effect of a decrease of signal is to 
increase the relative SM contribution in the samples. Overall, the values of the asymmetries 
will therefore be more SM-like, i.e., the central values will shift towards the value obtained 
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Figure 4: Total cross section in fb (sig- 
nal plus background) on peak \Mtj — 
M z >\ < 100 GeV, for M z , = 2.5 TeV 
for Z' — > e + e~, bb and ti, in the (g[, 
7j) plane, for decoupled neutrino (solid 
lines) and for m„ h = 50 GeV (dashed 
lines). The benchmark models are high- 
lighted as (red) dots on the (dotted) sce- 
nario lines. 
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in the SM. Obviously, this can happen only when the signal-to-background ratio is altered: 
in the case of leptons, where the SM contribution is negligible, the shift of the central values 
of the asymmetries is negligible. 

Comparing figure [I] to figure [2] and to figure [4], it is clear that the shape of the current 
exclusion limits is driven by the BR of the Z' boson into charged leptons: the limits are 
weaker where the BRs are smaller, i.e., in the region between the U(l) x and the ^ scenarios. 
It is also clear that searches performed in the top-quark final states do not improve this 
behaviour, being the BK(Z' — > ti) minimised near the U(l) x line. On the contrary, the 
6-quark final state's BR and cross sections are maximised here, so we checked whether 
exclusions derived in this final state could improve the overall results. It turns out that the 
sensitivity in this channel, even for comparable integrated luminosity, i.e., 5 fb -1 as in the 
analysis of Ref. [16J, is not comparable to the one exploiting the much cleaner di-lepton final 
state already discussed, that therefore yields the tighter constraints in the whole parameter 
space. 

3.4 Asymmetries 

We have so far described the features of the minimal Z' model in the (g[, g) plane by 
analysing standard variables, such as BRs, total width and event rates. For a complete 
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profile of the model, we now move on to study the asymmetries at the Z' peak. Performing 
a scan over the gauge couplings, the integrated values of asymmetries are presented in the 
(g[, g) plane for all three final states (Z 1 — > £ + £~ ,ti,bb). The observables were computed 



using the code described in sect. 2.4 implementing the cuts in eqs.(19)-(20) and folding 
in the relevant reconstruction efficiencies discussed in sect. 12.4.41 in order to determine the 
statistical uncertainties plotted underneath (always for 100 fb _1 of integrated luminosity). 
Finally, we distinguish among light leptons £ = e, /z, to be used in evaluating Arfb, and 
taus, essential to measure Al in the leptonic sector. 

We present here our results for Arfb and Al as representatives for charge and spin 
variables, respectively. Although we have produced results for it, we will not present A LL , 
because of its capabilities in the top final state only, while the aim of this paper is to compare 
predictions in different final states to assess the distinguishability of models. 

3.4.1 Charge asymmetry: Arfb 

Figure [5] shows the rapidity dependent forward-backward asymmetry along with its statistical 
uncertainty in the chosen three final states. In this case, the electron and muon final states 
can be used and possibly combined to reduce the statistical uncertainty further. One can 
clearly see the asymmetry vanishing in each final state along the scenario line corresponding 



to a zero value of one of the chiral couplings (vector or axial), as described in sect. 2.3 



The magnitude of the statistical uncertainties matches the total cross section plots in 



figure |4j in accordance with eq. (30) and, outside of the cases near the trajectories where a 
particular final state has vanishing asymmetries, the uncertainties can be as low as around 
25% in the case of the top final state. In the leptonic final state, the uncertainties are 
comparatively smaller, apart form the large spike at the SM point due to the lack of cross 
section. Finally, the bb final state still appears to perform reasonably well with errors of 
order 30-50%, considering the extremely low reconstruction efficiency. The sensitivity to 
the relative sign of the couplings in the final state is blurred due to the dependence on the 



product of the initial and final state couplings, as discussed in sect. 2.4.1 

The overall assessment of the visibility of the Z' boson suffers from the large uncertainties 
(especially if compared to i^, as we will see later). As shown in figure [6j where also the 
coloured projection labels the final state which offers the highest significance for each pair 
of coupling values is given, 100 fb _1 are enough to gain a sensitivity equal or greater than 3 
almost everywhere. The dominance of the leptonic final state is remarkable, this final state 
being sufficient to cover most of the parameter space. This is due to the large positive value 
of this asymmetry for the SM, while in our model the leptonic final state yields a value for 
Arfb much smaller or even slightly negative in most of the parameter space. It is peculiar 
that, for large negative values of g, where also A e RFB is large, the already small errors are not 
sufficiently small to allow one to distinguish the Z' boson from the SM. Nonetheless, a good 
discrimination power is there provided by top quark final states, although with significance 
always below 5 due to the larger uncertainties. 

The effect of introducing the heavy neutrinos is to reduce significances across the board, 
which results in a slight enlargement of the grey area in which a 3a significance cannot be 
obtained in any final state. The reason for this is twofold, although driven by a single cause: 
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Figure 5: Arfb on peak \M^j — 
M z >\ < 100 GeV, for M z , = 2.5 
TeV for Z' — > e + e~, bb and ti in 
the (g[, g) plane. Statistical uncer- 
tainties are shown underneath each 
plot in the same plane assuming 100 
fb _1 of integrated luminosity. The 
benchmark models are highlighted as 
(green) dots on the (dotted) scenario 
lines and, where relevant, a contour 
shows where the asymmetry vanishes. 



the reduction of the signal rates in the considered final state. First, as discussed at the end 



of section 3J3 the central value of the asymmetry is shifted towards SM values due to a 
reduced signal-over- SM ratio of events. This does not affect the lepton final state, which 
significance is reduced due to the second effect: smaller total rates mean larger errors, which 
naively reduce the significance. 

3.4.2 Polarisation asymmetry: Al 

Figure [7] shows the single polarisation asymmetry along with its statistical uncertainty in 
the usual three final states. Recall that in the case of the leptonic final state labelled ee, 
the measure would be obtained from the tau final state with the discussed reconstruction 
efficiency times BR of 10%. Once again and even more clearly than for Arfb, one can 
see the asymmetry vanishing in each final state where expected. The change in sign of 
the asymmetry at different points in parameter space corresponds to a change in relative 
sign of the vector and axial couplings of the Z' boson to the final state as discussed in 



sect. 2.4.2 This is most pronounced in the di-lepton case because of the comparatively 
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Figure 6: Significance of A RFB with respect to the SM prediction on peak \Mjj— M z > \ < 100 
GeV, for Mz< = 2.5 TeV for Z' — > e + e~ (red), bb (blue) and ti (green) in the (g[, g) plane 
at the LHC at 14 TeV both for decoupled heavy neutrinos (left) and for heavy neutrinos 
of 50 GeV (right). The coloured projection denotes the final state which offers the most 
significance for each pair of coupling values. We assume here a luminosity of 100 fb _1 . 

low SM background, which dilutes the asymmetry in the other final states. The behaviours 
of the statistical uncertainties do not greatly differ from those in the Arfb plots as they 
are largely determined by the cross section values, although they are about a factor of two 
smaller, making this variable the one with the greatest discrimination power. 

Overall this gives an indication that all channels could be useful in probing the full 
parameter space of such a Z' model in complementing one another by having different areas 
where the sensitivity is the best. This is shown in figure [8] where significances in excess 
of 3a appear over almost all of the (g[, g) plane. Although leptons almost always provide 
the best significance, due to the large negative SM value as compared to the large positive 
signal values over the parameter space, in the areas where the two are similar for the lepton 
final state, bb and ti come into play and provide adequate 'coverage'. Again, having coupled 
heavy neutrinos results in an overall reduction of the signal significance with respect to the 
SM and in a slight enlargement of the grey area in which 3a significance cannot be obtained 
in any final state. However, comparing figure [6] to figure [8j it is clear how these grey areas 
overlap only in the region where both gauge couplings g[ and g are small, showing a great 
complementarity of observables and final states. 

So far we have discussed how well the asymmetry produced by a Z' boson in our model 
can be distinguished from the SM background. One would also like to disentangle different 
combinations of gauge couplings, in turn leading to their absolute measurement. Although 
the full analysis and especially this last part are beyond the scope of this paper, we would 
like to briefly discuss the first point. 

Figure |9] shows that all our benchmark models can well be discriminated with 100 fb _1 
of data. Especially, degenerate models in one final state (such as B — L and x m the case of 
top quarks) are well separated either using leptons or, eventually, 6-quarks. The last frame, 
displaying the t- versus 6-quark case, clearly shows the impact of including heavy neutrinos: 
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Figure 7: Al on peak \Mfj — Mz>\ < 
100 GeV, for M z > = 2.5 TeV for Z' -> 
e + e~, 66 and tt in the g[, 'g plane at 
the LHC at 14 TeV. Statistical uncer- 
tainties are shown underneath each 
plot in the same plane assuming 100 
fb _1 of integrated luminosity. The 
benchmark models are highlighted as 
(green) dots on the (dotted) scenario 
lines and, where relevant, a contour 
shows where the asymmetry vanishes. 

represented by dashed crosses, the asymmetries with light are closer to the SM values and 
have larger errors. 
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Table 3: Significance for A L for the LHC at 14 TeV for 100 fb" 1 and M z > = 2.5 TeV for 
the common benchmark points in the r, t, b final states. Upper triangle for decoupled heavy 
neutrinos and lower triangle for m^ h =50 GeV. 

In table [3] we collect the significance of each benchmark point with respect to any other 
when Al is measured in each of the final states, would the measure in 6-quark final state 
be available, both for decoupled neutrinos (upper triangle) and for m Ufi = 50 GeV (lower 
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Figure 8: Significance of A L withe respect to the SM prediction on peak \Mfj — M Z '\ < 100 
GeV, for Mz> = 2.5 TeV for Z' — > e + e~ (red), bb (blue) and ti (green) in the (g[, g) plane 
at the LHC at 14 TeV both for decoupled heavy neutrinos (left) and for heavy neutrinos 
of 50 GeV (right). The coloured projection denotes the final state which offers the most 
significance for each pair of coupling values. We assume here a luminosity of 100 fb -1 . 
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Figure 9: On-peak values with 
relative errors for Mz> = 2.5 TeV in 
different combinations of final states, 
for the selected benchmark models in 
sect. I2.3I Solid crosses are for decou- 
pled neutrinos, while dashed crosses 
are for m Vh = 50 GeV. Leptons do 
not show any appreciable difference, 
hence they are not displayed. 
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triangle). We first notice that most of the models are quite distinguished from each other 
just by looking at taus, that generally speaking delivers higher discrimination power due to 
the smaller errors. The only exception is the separation between fu and B — L, both with 
vanishing asymmetries in the lepton final state. Here, the supplementary information from 
the top final state already proves to be sufficient to discriminate among all the models. If the 
inclusion of tops already helps to disentangle the models, it is clear that were it available, 
the measure in 6-quarks would be of great help to fully distinguish them: the discrimination 
power is always above the 4<r level in at least one final state, the 6-quarks being especially 
relevant to discriminate between fj and p. 

Regarding the inclusion of heavy neutrinos, their impact is that of reducing the absolute 
value of the asymmetry and to increase its statistical error. Taus are not affected by their 
presence and therefore show the same significance. As well, the absolute values of the 
asymmetries for B — L and p are not affected by heavy neutrinos, given that the former 
has vanishing asymmetries everywhere and that the latter has negligible BRs into them. 
Therefore, their disentanglement from any other model can only get worse, due to the reduced 
central values of the latter. The same is true also for x when tops are considered, given that 
its asymmetries in this final state are vanishing. When instead models with finite values of 
the polarisation asymmetry in a specific final state are compared (i.e., R and using tops 
and R, fj and x using 6's), the discrimination power therein gets enhanced when considering 
coupled heavy neutrinos. This seems counterintuitive, but it can be understood looking at 
figure [2]in terms of BRs. In fact, the central values change differently and the various models 
considered here are better separated, despite the slightly larger statistical errors. 

As a last example, we want to make the case of whether all possible gauge coupling 
combinations could be resolved in our approach. It is not hard to see in figure [8] that there 
exist combinations of couplings that yield similar (within errors) asymmetries in top and tau 
final states. Possibly, the study of differential and angular distributions of charged leptons, 
as in Ref. [T5], can fully resolve the whole parameter space. However, here we would like 
to make a different point: if the measure of the (spin) asymmetries in 5-quark final states 
becomes available, it would be of great interest. To this purpose, we select four pairs of 
gauge couplings, shown in figure [lO^top-left) that have degenerate asymmetries both for 



leptons and for tops (see fig. 10 [top-right)). The pairs of points labelled "1" and "2" are 
examples of combinations that could be discriminated when the asymmetry in 6-quarks is 
also considered. However, this is still not possible over the whole parameter space, as clear 
from the points labelled "3" and "4" , fully degenerate, so in any case one would require to 
integrate such information with a study of standard observables such as (differential) cross 
sections, lineshape fits and BRs, reinforcing the importance of making use of all accessible 
information. We remark that these last statements are valid at 100 fb^ 1 and for the efficiency 
factors employed. Higher amounts of data, and better efficiencies, will certainly improve the 
discrimination power, as well as new generations of micro-vertex detectors could make the 
6-quark final state competitive with the r-lepton and t-quark ones. 
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Figure 10: (Top-left) combinations of gauge couplings with some degenerate asymmetries. 
(Other frames) On-peak Al values with relative errors for combinations of final states for 
M z > =2.5 TeV and decoupled heavy neutrinos. 

4 Conclusions 

In summary, we have proven the feasibility of profiling a Z' boson, possibly discovered at the 
14 TeV LHC, thanks to the exploitation of the decays of the new gauge boson into t + t~, 
ti and, possibly, hh states, with respect to the time-honoured studies of DY channels only 
(i.e., decays into e + e~ and fi + fi~ final states). In fact, the former signatures afford one with 
the possibility of not only defining standard charge asymmetries, the only ones accessible 
in the latter, but also spin ones. In essence, such heavy fermions (i.e., r, b and t states) 
can transmit their spin properties to their decay products rather efficiently (even though 
just in some specific cases), thereby making it possible to define single- and/or double-spin 
asymmetries that are peculiarly dependent on the quantum numbers of the new gauge boson. 
Their measure is certainly feasible for tau and top-quark final states, while, at this moment 
in time, the 6-quark final state may prove difficult, if not possible at all, to analyse due 
to experimental limitations in reconstruction, as well as to the highly boosted kinematics. 
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Further, based on a dedicated parton level simulation, yet including some selection criteria 
and realistic detector efficiencies, we have argued that in our approach neither DY nor any of 
the above new channels can be used alone to fully probe all the parameter space. Rather, one 
way to render this possible is to combine two or more of these channels, including also spin 
observables, as the aforementioned final states show different sensitivity to the observables 
studied. Nonetheless, due to the continuous nature of the parameter space, degenerate 
regions still exist. We finally made also the point on how the presence of light extra states 
(the heavy neutrinos) can alter the observables under study, and that such impact could be 
resolved. Such conclusions can be extended to a large variety of the most popular Z' models 
studied nowadays, all compliant with current experimental limits, and in presence of SM 
irreducible background effects. 

Ultimately, the goal of profiling an observed resonance would be to measure all the 
parameters of its fermionic interactions. In the minimal CcLSG, clS discussed in section 2, there 
are 5 independent couplings. Following Ref. [15], at least one more linearly independent 
observable is required than those available in the light lepton sector, where some degeneracy 
still ocours. Considering multiple final states and their spin asymmetries could provide 
the necessary measurements to extract the fermionic couplings in this minimal case, that 
become essential when testing the univerality of couplings or when considering models with 
more general coupling structures. Finally, our continuous scanning of the parameter space 
highlighted regions where leptonic final states are not the most sensitive ones, hence where 
alternative final states would prove to be more effective. 

Acknowledgements 

We are grateful to A. Belyaev for the suggestion of the pr{b) cut and to I.R. Tomalin for 
pointing to us the LEP measurements of the polarisation of b-hadrons. LB thanks Stan Lai 
for helpful discussions about taus. Finally we would like to thank the anonymous referee 
for referring us to [IE] and additionally stimulating a great deal of internal discussion and 
development of our understanding of polarisation asymmetries and their role in investigating 
the Z' couplings. The work of KM and SM is partially supported through the NExT Institute. 
LB is supported by the Deutsche Forschungsgemeinschaft through the Research Training 
Group grant GRK1102 Physics of Hadron Accelerators. 

References 

[1] For a review, see M. Schmaltz and D. Tucker-Smith, Ann. Rev. Nucl. Part. Sci. 55, 
229 (2005). 

[2] For a review, see J. L. Hewett and M. Spiropulu, Ann. Rev. Nucl. Part. Sci. 52, 397 
(2002). 

[3] For a review, see: R. N. Mohapatra, Unification and Supersymmetry, Springer, New 
York, 1986. 



24 



[4] For a review, see J. L. Hewett and T. G. Rizzo, Phys. Rept. 183, 193 (1989). 

[5] For a review, see D. J. H. Chung, L. L. Everett, G. L. Kane, S. F. King, J. D. Lykken 
and L. T. Wang, Phys. Rept. 407, 1 (2005). 

[6] M. J. Strassler and K. M. Zurek, Phys. Lett. B 651, 374 (2007). 

[7] CDF Collaboration, Phys. Rev. Lett. 102, 091805 (2009) and Phys. Rev. Lett. 102, 
031801 (2009); DO Collaboration, note 4375-Conf and note 4577-Conf. 

[8] C MS Collabora tion, Phys. Lett. B 714, 158 (2012); ATLAS Collaboration, 
|arXiv:1209.2"535] [hep-ex] ; 

[9] A. Papaefstathiou and K. Sakurai, JHEP 1206, 069 (2012). 

[10] L. Basso, K. Mimasu and S. Moretti, JHEP 1209, 024 (2012). 

[11] D. Choudhury, R. M. Godbole, S. D. Rindani and P. Saha, Phys. Rev. D 84, 014023 
(2011); R. M. Godbole, K. Rao, S. D. Rindani and R. K. Singh, JHEP 1011, 144 

(2010) ; B. Xiao, Y. K. Wang, Z. Q. Zhou and S. H. Zhu, Phys. Rev. D 83, 057503 

(2011) and Phys. Rev. D 83, 094022 (2011); E. L. Berger, Q. H. Cao, C. R. Chen and 
H. Zhang, Phys. Rev. D 83, 114026 (2011); R. Diener, S. Godfrey and T. A. W. Martin, 
Phys. Rev. D 83, 115008 (2011). 

[12] T. Stelzer and S. Willenbrock, Phys. Lett. B 374, 169 (1996); G. Mahlon and S.J. 
Parke, Phys. Rev. D 53, 4886 (1996) and ibidem 81, 074024 (2010); W. Bernreuther, 



A. Brandenburg, Z.G. Si and P. Uwer, Phys. Rev. Lett. 87, 242002 (2001) and hep- 



ph/0410197 



[13] See, e.g.: S. Jadach and Z. Was, in L Z Physics at LEP1', edited by G. Altarelli, R. 
Kleiss and C. Verzegnassi, preprint CERN 89-08. 

[14] ATLAS Collaboration, Eur. Phys. J. C 72, 2062 (2012); ATLAS-CONF-2011-152. 

[15] F. Petriello and S. Quackenbush, Phys. Rev. D 77 (2008) 115004. 

[16] CMS-PAS-EXO-11-008. 

[17] J. Qlzem [CMS Collaboration], JINST 6, C12039 (2011), 

http : //indico . cern. ch/ contributionDisplay .py?contribId=39fesessionId=9& 
confld=190941. 

[18] B. Mele and G. Altarelli, Phys. Lett. B 299, 345 (1993); G. Bonvicini and L. Randall, 
Phys. Rev. Lett. 73, 392 (1994); A. F. Falk and M. E. Peskin, Phys. Rev. D 49, 3320 
(1994); M.M. Tung, J. Bernabeu and J. Penarrocha, Phys. Lett. B 418, 181 (1998). 

[19] D. Buskulic et al. [ALEPH Collaboration], Phys. Lett. B 365, 437 (1996); G. Abbiendi 
et al. [OPAL Collaboration], Phys. Lett. B 444, 539 (1998); P. Abreu et al. [DELPHI 
Collaboration], Phys. Lett. B 474, 205 (2000). 



25 



[20] E. Accomando, A. Belyaev, L. Fedeli, S. F. King and C. Shepherd-Themistocleous, 
Phys. Rev. D 83, 075012 (2011). 

L. Basso, S. Moretti and G. M. Pruna, JHEP 1108, 122 (2011). 
F. del Aguila, M. Masip and M. Perez- Victoria, Nucl. Phys. B 456, 531 (1995). 
M. Carena, A. Daleo, B.A. Dobrescu and T.M.P. Tait, Phys. Rev. D 70, 093009 (2004). 
P. H. Chankowski, S. Pokorski and J. Wagner, Eur. Phys. J. C 47, 187 (2006). 
A. Ferroglia, A. Lorca and J. J. van der Bij, Annalen Phys. 16, 563 (2007). 

E. Salvioni, G. Villadoro and F. Zwirner, JHEP 0911, 068 (2009). 

F. del Aguila, G. D. Coughlan and M. Quiros, Nucl. Phys. B 307, 633 (1988) [Erratum- 
ibid. B 312, 751 (1989)]. 

T. Appelquist, B. A. Dobrescu and A. R. Hopper, Phys. Rev. D 68, 035012 (2003). 

P. Abreu et al. [DELPHI Collaboration], Z. Phys. C 65, 603 (1995). 

H. Murayama, I. Watanabe and K. Hagiwara, KEK Report 91—11, January 1992. 

T. Stelzer and W. F. Long, Comput. Phys. Commun. 81, 357 (1994). 

J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 
07, 012 (2002). 

G. P. Lepage, J. Comp. Phys. 27, 192 (1978) [Erratum, preprint CLNS-80/447, March 
1980]. 

A. Pukhov, arXiv:hep-ph/0412191| 
A.V. Semenov, j arXiv:hep-ph/9608488| 

N. D. Christensen and C. Duhr, Comput. Phys. Commun. 180, 1614 (2009) 1614; 
N. D. Christensen, P. de Aquino, C. Degrande, C. Duhr, B. Fuks, M. Herquet, F. Mal- 
toni and S. Schumann, Eur. Phys. J. C 71, 1541 (2011). 



M. Bondarenko et at, in arXiv: 1203. 1488 [hep-ph]; see also 
https://hepmdb.soton.ac. uk/. 

See https : //f eynrules . irmp .ucl . ac .be/wiki/B-L-SM. 

R. W. Brown, D. Sahdev and K. O. Mikaelian, Phys. Rev. Lett. 43, 1069 (1979); 
J. H. Kuhn and G. Rodrigo, Phys. Rev. D 59, 054017 (1999). 

[40] M. Arai, N. Okada, K. Smolek and V. Simak, Acta Phys. Polon. B 40, 93 (2009). 



26 



[41 
[42 

[43 
[44; 

[45 
[46 
[47 



K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B 274, 1 (1986). 

J. L. Hewett, J. Shelton, M. Spannowsky, T. M. P. Tait and M. Takeuchi, Phys. Rev. 
D 84, 054005 (2011). 

E. Alvarez, Phys. Rev. D 85, 094026 (2012). 

CMS Collaboration, CERN-LHCC-2006-001; ATLAS Collaboration, CERN-LHCC- 
99-14 and CERN-LHCC-99-15. 

M. Frank, A. Hayreter and I. Turan, Phys. Rev. D 84, 114007 (2011). 
CMS PAS BTV-11-004. 

See, e.g.: https : //twiki . cern. ch/twiki/bin/view/CMSPublic/SWGuideCategory 
BasedElectronID#Electron_Ef f iciency . 



[48] L. Basso, O. Fischer and J. J. van der Bij, arXiv: 1207.3250 [hep-ph]. 



27 



